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, We study orthogonal and symplectic matrix models with polynomial potentials and 

multi interval supports of the equilibrium measure. For these models we find the bounds 
(similar to the case of hermitian matrix models) for the rate of convergence of linear 
eigenvalue statistics and for the variance of linear eigenvalue statistics and find the log- 
arithms of partition functions up to the order O(l). We prove also universality of local 
eigenvalue statistics in the bulk. 

C3 ■ 

1 Introduction and main results 

In this paper we consider ensembles of random matrices, whose joint eigenvalues distribution 
is 

m 

J2: PnA^~An) = Q^iv]iie-^y* n \h-x J f = Q-y H ^-^ )/ \ a.i) 

CN \ i=l l<i<j<n 

■ where the function H, which we call Hamiltonian to stress the analogy with statistical me- 

chanics, and the normalizing constant Q n)J g[V] have the form 



H{Xi, ...,\ n ) = -nJ2 v ( x i) + Yl lo S l A ' - A il> 

i=l v£j 

QnAV] = [ e^- x ^ 2 dX 1 . . . dX n . (1.2) 



The function V, called the potential, is a real valued Holder function satisfying the condition 

F(A)>2(l + e)log(l + |A|). (1.3) 

This distribution can be considered for any (5 > 0, but the cases f3 = 1,2,4 are especially 
important, since they correspond to real symmetric, hermitian, and symplectic matrix models 
respectively. 

We will consider also the marginal densities of (jl.ip (correlation functions) 

Pi n }(X 1 ,...,\i) = Pn,/3(Ai, ...Xi,Xi + i, X n )d\i + i...d\ n , (1.4) 

Jm.™- 1 

and denote 

%{(•••)}= /(• • • K/3(Ai, X n )dX u . . . dX n . (1.5) 
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It is known (see [21 [12]) that if V' is a Holder function, then the first marginal density pf^ 
converges weakly to the function p (equilibrium density) with a compact support a. The 
density p maximizes the functional, defined on the class Ai\ of positive unit measures on R 

£ v (p) = max {L[dm,dm} - [ V(X)m(dX)\ = £[V], (1.6) 

meMi { J J 

where we denote 

L[dm,dm] = J log |A — p\dm(X)dm(p), L[f, g] = J log | A — p\f(X)g(p)dXdp. (1-7) 
The support a and the density p are uniquely defined by the conditions: 



as 



v(X) := 2 y log |/i - X\p(p)dp - V(X) = supv(X) := w*, X £ a 
v(X) < supv(X), X a, a = suppjp}. 

For (3 = 2 it is well known (see f!4j ) that all correlation functions (jl.4j) can be represented 
pg(A l5 ...,A,) = ^^&et{K n>2 {X j: X k )} l j>k=l , (1.9) 



where 

n-l 

K n>2 (X,p) = ^ n) (X)^ n) (p). (1.10) 

This function is known as a reproducing kernel of the orthonormalized system 

4 n) (X) = eM-nV(X)/2}p { l n) (X), I = 0,..., (1.11) 

in which {p[ n ^}f = Q are orthogonal polynomials on R associated with the weight w n (X) = 
e- nV ^,i.e., 

"p ( i n) (X)p^(X)w n (X)dX = 6 l>m . (1.12) 

The orthogonal polynomial machinery, in particular, the Christoffel-Darboux formula and 
Christoffel function simplify considerably the studies of marginal densities (|1.4|) . This allows 
to study the local eigenvalue statistics in many different cases: bulk of the spectrum, edges 
of the spectrum, special points, etc. (see [IB], [17], [5], [5], [3], [15], [11]). 

For j3 = 1,4 the situation is more complicated. It was shown in [22] that all correlation 
functions can be expressed in terms of some matrix kernels (see (|1.13p - f|l . lTj) below). But 
the representation is less convenient than (|1 .9j) - (jl.lOp . It makes difficult the problems, 
which for (3 = 2 are just simple exercises. For example, the bound for the variance of linear 
eigenvalue statistics (|1.2ip for /3 = 1,4 till now was known only for one interval a (see [12]). 
while for /3 = 2 it is a trivial corollary of the Christoffel-Darboux formula for any a. 

The matrix kernels for /3 = 1,4 have the form 

*-'<»■"> - U^^a - ,» ~iXf) for " = eve - (L13) 
^-H^tii) (L14) 
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where 

n-l 

5 nil (A,/x) = - £ ^(A)^)):^^)^), (1.15) 

71-1 

5 nA 4(A, M ) = - £ (^ n) )'(A)Pi n) )- fc Vi n) (/x), (1.16) 

i,fc=0 

e(A) = ^sgn(A), sgn denotes the standard signum function, 

(e/)(A) := / e(A -//)/(//) d/idA'. 

Z?i n) and i\4 n) in (fTT5|) and (fTT6j) are the left top corner n x n blocks of the semi-infinite 
matrices that correspond to the differentiation operator and to some integration operator 
respectively. 

^ ) :=((^ ) )',V'i n) ) , D^ = {D^}l k=Q , (1.17) 

Both matrices D<£? and AT^ are skew-symmetric, and since e(^>j )' = > we have for any 
j, I > that 

oo 

8 3l = {e{i,f ) )\^ l ) = Y J {Dt ) ) 3 k{M^) kl DL n) MW = l = MWDW. 

fc=0 

It was observed in [23] that, if V is a rational function, in particular, a polynomial of degree 
2m, then the kernels S n x, Sn4 can be written as 

2m- 1 

S njl (A, M ) = ^n,2(A,/x)+n £ F^^U^nlk^ ( L18 ) 

j,fc=-(2m-l) 
2m-l 

S n/2A (\,v) = K n , 2 (\,n) + n £ ^ 4 Vi"i-(^)^S fc (/x), 

i,fc=-(2m-l) 

where i^^j ^ij! can be expressed in terms of the matrix T" 1 , where T n is the (2m — 1) x 
(2m — 1) block in the bottom right corner of Mn , i-e., 

(T n ) jk := (L»(")M^) n _ 2m+iin _ 2m+fe , 1 < j, fc < 2m - 1. (1.19) 

The main technical obstacle to study the kernels S nt i,S nt 4 is the problem to prove that 
(T~ W are bounded uniformly in n. Till now this technical problem was solved only in a 
few cases. In the papers [HIE] the case V(X) = A 2m (l + o(l)) (in our notations) was studied 
and the problem of invertibility of T n was solved by computing the entries of T n explicitly. 
Similar method was used in [8] to prove bulk and edge universality (including the case of 
hard edge) for the Laguerre type ensembles with monomial V . In the paper [21] the problem 
of invertibility of T n was solved also by computing the entries of T n for V being an even 
quatric polynomial. In |19[ [20] similar problem was solved without explicit computation of 
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the entries of T n , it was shown that for any real analytic V with one interval support of the 
equilibrium density (M,l ) —1 is uniformly bounded in the operator norm. This allowed us to 
prove bulk and edge universality for j3 = 1 in the one interval case. 

But there is also a possibility to prove that T n is invertible with another technique. As a 
by product of the calculation in [22] one also obtains relations between the partition functions 

Q n fi and the determinants of Mn and '■ 

where T n := Ilj=o > ano - 1j ^ s ^ ne lading coefficient of pj 1 \\) of (| 1 . 1 2 1) . It is also known 
(see [H]) that Q nj2 = T 2 Jn\. Since d£ ] m£ } = 1 and (D&) jk = for \j - k\ > 2m - 1 (see 
(pT5]) ). we have = 1 + A n with A n being zero except for the bottom 2m — 1 rows, 

and we arrive at a formula, first observed in [21] : 

de t (T„,^et( D .».A4"',^(|^) 2 . (1.20) 

Hence to control det(T n ), it suffices to control log Q n> p for /3 = 1, 2, 4 up to the order 0(1). In 
the paper [15] the corresponding expansion of log Q n ,/3 was constructed by using some gener- 
alization of the method of [12]. The original method was proposed to study the fluctuations 
of linear eigenvalue statistics 

n 

AfnM = X>(Ai), (1.21) 
i=l 

in particular, to control the expectation and the variance of n~ l N n [tp\ up to the terms 0(n~ 2 ) 
for any /?, but only in the case of one interval support of the equilibrium measure and 
polynomial V, satisfying some additional assumption. The method was used also to prove 
CLT for fluctuations of M n [(f). In the paper [15] the method of p2] was simplified, that 
allowed to generalize it on the case of real analytic V with one interval support of p, without 
any other assumption. Unfortunately, there is no hope to generalize the method of [121 [T5] 
on the case of multi interval support a directly, because the method is based on solving of 
some integral equation (see Eq. (|2.12p below) which is not uniquely solvable in the case of 
multi interval support. 

In the present paper the problem to control Qn./sIX] for f3 = 1,2,4 is solved in a little 
bit different way. We prove that for analytical potential V with g-interval support a of the 
equilibrium density Q n ^ can be factorized to a product of Qk^^Wa^}, a = 1, . . . ,q, where 

k* a ~ /x*n (see (|1.33j) ). and the "effective potentials" Va°^ (see (|1.35j) ) are defined in terms of 
cr, V and p. 

To be more precise let us formulate our main conditions. 

Condition CI. V is a polynomial of degree 2m with positive leading coefficient, and the 
support of its equilibrium measure is 

q 

o- = (J a a , a a = [E 2a ^,E 2a } (1.22) 

a=l 

Condition C2. The equilibrium density p can be represented in the form 

P(A) = 7 ^-P(A)9X 1 / 2 (A + iO), inf |P(A)| > 0, (1.23) 

Z7T Aecr 
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where 

2q 

X(z)=H(z-E a ), (1.24) 

a=l 

and we choose a branch of X 1 ^ 2 (z) such that X l l 2 (z) ~ z q , as z — >■ +00. Moreover, the 
function v defined by [1.8]) attains its maximum only if A belongs to a. 



Remark 1 It is known (see, e.g., fljj) that for analytic V the equilibrium density p has the 
form [1.23]) - [1.2$ with P > 0. The function P in [1.23]) is analytic and can be represented 
in the form 

Hence condition C2 means that p has no zeros in the internal points of a and behaves like 
square root near the edge points. This behavior of V is usually called generic. 

We will use also the notations 

q 

a z = U <Ta > £ ' cr ». e = t^a-l ~ E 2a + e], (1.26) 
0=1 

dist {a a , e , o- a ',e} > 5 > 0, a ^ a . 

The first result of the paper is the theorem which allows us to control log Q n ^ in the one 
interval case up to 0(1) terms. Since the paper [2] it is known that 

logQ^y] = f3n 2 £[V}/2 + 0(n log n). 

But, as it was discussed above, for many problems it is important to control the next terms 
of asymptotic expansion of log Q n ,/3 (see also the discussion in [jJJ, where the expansion in re -1 
was constructed for (3 = 2 and V being a polynomial, close in a certain sense to Vo(A) = A 2 /2.) 

We would like to note that almost all assertions of Theorem [T] below were obtained in 
|15j . The difference is that here we need to control that the remainder bounds are uniform in 
some parameter 77, which we put in front of V in H (see (jl.2p ). Note also that it is important 
for us that here V may be non polynomial function analytic only in a some open domain 
D C C containing a. 

We will use below the Stieltjes transform, defined for an integrable function p as 

»M = / 4^- (1-27) 



Theorem 1 Let V satisfy [1.3]) . the equilibrium density p (see [1.8]) ) have the form [1.23]) 
with q = 1, and a = supp/j = [a, b]. Assume also that V is analytic in the domain D D u £ . 
Consider the distribution [1.1]) with V replaced by 77V. Then there exists e\ > such that for 
any rj : \r] — 1| < e\ we have: 

(i) The Stieltjes transform gn (z) [1.27\) of the first marginal p[j^ of [1.1]) for z such that 
d(z) := dist{z,o" e } > n -1 / 6 logn has the form 

g^(z)=g v (z)+n~ 1 u n>v (z), (1.28) 
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where g v (z) is the Stieltjes transform of the equilibrium density p v , which maximizes 8 [rjV] 
of il.6\) , X v of [1.2$ corresponds to the support [a^ , 6^] of p v , and P^ is defined by $1.25]) 
for rjV . The contour C here is chosen sufficiently close to [a^,^] to have z and all zeros of 
P v outside of it. The remainder bound is uniform in \rj — 1| < e\. Moreover, for any ip with 
bounded fifth derivative 



l p(X)(p^f(X)-p v (X))dX = n- 1 0(\\ l p\\ 0O + ||^ 6) ||oo), (1-29) 



Halloo := sup \<p(\)\. 



(ii) There exists an analytic in~D\a e function u* such that 

u n , v (z)-u n , 1 (z) = (r l -l)u*{z) + 0((r ] -l) 2 ) + 0(n- 1 ), \$fz\ >d. (1.30) 
(Hi) If Q^ a is defined by hl.2\) for rjV , then 



l °&Qn'/3 = l °gQn,p + — £W\ + — g~ +n\l - -J log- 



'/ 



dt d> ( V V(z)-V v ^(z))u n , v (z,t)dz, 



_[3n 
2 

(2//3-1) r g' v ((,t)d( 

u n ,r,(z,t) = f ( — + 0(n 1 ), 1.31 

where £[r]V], X„ and [ou,6J are the same as in (ii), 

V^\z)=2{z-c v f/d 2 r} , c v = (a v + b v )/2, d v = {\ - a v )/2, 
P V (X, t) = tP v (X) + fc^, g v (z, t) = tg v (z) + ^~w^(z - c„ - 

Q* n p is defined by the Selberg formula 



(f) ■ (L32) 

and the remainder bound is uniform in \rj — 1| < ex- 

The next theorem establishes some important properties of the symplectic and orthogonal 
matrix models, in particular, it gives the bound for the rate of convergence of linear eigenvalue 
statistics and the bound for their variances. In order to formulate the theorem, we define for 
even n 

H* a = [ p(X)d\, k* a := [np* a ] + d a , (1.33) 



where [x] means an integer part of x, and we have chosen d a = 0, ±1,±2 in a way which 
makes fc* even and 

For each a a ^ £ we introduce the "effective potential" 

V^ a) (X) = l CTa , e (A)(y(A)-2 / \og\X -p\p(p)dp), (1.35) 
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and denote X* the "cross energy" 

£*: = £ / dA / d//log|A-M|p(A)p(M). (1-36) 

/ „ / J (Try »/ <7„,/ 



Theorem 2 If potential V satisfies conditions C1-C2 and n is even, then the matrices 
and in (TIM) are bounded in the operator norm uniformly in n. Moreover, for any 
smooth ip and f3 = 1,2,4 we have 



^(A)(pg(A)-p(A))dA 



< -Halloo, (1-37) 
n 



E 



MnVp] ~ Ep{N n [<p]} 2 \ < C\\ip' u2 



where ||.||oo is defined in M.29j) . The logarithm of the normalization constant Qn^W] can be 
obtained up to 0(1) term from the representation 

log(Q n AV}/n\) = X)log(Q fc ^[^ a )]/fc*!)-^E* + 0(l), (1.38) 

a=l Z 

where Vd and X* are defined in hi. 3 5]) and hl.36\) . 

As it was mentioned above, Theorem [2] together with some asymptotic results for orthogonal 
polynomials of [5] may be used to prove universality of local eigenvalue statistics of matrix 
models (jl.lj) . In order to state our theorem on bulk universality, we need some more notations. 
Define 

T ^ , . sin7rt 
Koo{t) :-- 




irt ' 

jy Koo (t)dt-e(Z-r)) K^ti-OJ' 

J* "Kooftdt Koch-®)' 

Furthermore, we denote for a 2 x 2 matrix A and A > 

Theorem 3 Let V satisfy conditions C1-C2. Then we have for (even) n — > oo, Xq 6 E with 
p(X ) > 0, and for (3 G {1,4} that 

q- l Knf (Ao + H/Qn, A + v/q n ) = K<g>(S, rj) + 0(n~ 1 / 2 ), 
q- l K^ A (A + t/q n , A + Vhu) = K^{i, rj) + 0(n-^ 2 ), 

where q n = np(Xo). The error bound is uniform for bounded £, rj and for Xq contained in 
some compact subset of U^ =1 (-E^a-i ; £?2a)- 

It is an immediate consequence of Theorem [3] that the corresponding rescaled /-point 
correlation functions 

P$ ( A o + £i/q n , ■ ■ ■ , Ao + ii/q n ) , p|,4 /2) (A + £i/q n , • • • , A + ii/q n ) 

converge for n (even) — > oo to some limit that depends on (3 but not on the choice of V. 

The paper is organized as follows. In Section [2] we prove Theorem [TJ In Section [3] we 
prove Theorems [2] and [3] modulo some bounds which are obtained in Section [H 
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2 Proof of Theorem [T] 



Proof of Theorem [7J 

(i). Take n-independent e, small enough to provide that a £ C D. It is known (see e.g. [2]) 
that if we replace in (jl.ip . (|1.4p . and (jl.5p the integration over R by the integration o~ e , then 
the new marginal densities will differ from the initial ones by the terms 0(e~ nc ) with some c, 
depending on e, but independent of n. Since for our purposes it is more convenient to consider 
the integration with respect to o~ e , we assume from this moment that this replacement is made, 
so everywhere below the integration without limits means the integration over a £ . 

Following the idea of |12| . we will study a little bit modified form of the joint eigenvalue 
distribution, than in (jl.ip . Namely, consider some real smooth function h(X) and denote 

V h>v (C) = vV(0 + -HC). (2.1) 

n 

Let Vn,fi,h-, ^/3,h{- ■ ■ }; P^Bh be the distribution density, the expectation, and the marginal 
densities defined by (11. ip . (11.5p . and (II. 4|) with V replaced by V^. 
By (jl.ip the first marginal density can be represented in the form 

/n 
e-^(A)/2TJ| A _ A .|^ e -^y M (A,)/2 Yl \X i -X j fd\ 2 ...d\ n . 
i=2 2<i<j<n 

(2.2) 

Using the representation and integrating by parts, we obtain 

z — X pn J (z—Xy n J (z — X){X — fi) 

Here 0(e~ nc ) is the contribution of the integrated term. In fact all equations below should 
contain 0(e _nc ), but in order to simplify formula below we omit it. 

Since the function P^pii^, /-0 is symmetric with respect to A, fi, we have 

p£pM X > M)dA<i/x _ f p^pl(X, fi)dXd^ f p^p^iX, fi)dXdfi _ f P^ h {X, fi)dXdfi 



(z-X)(X-fi) J ( z -X)(X-fi) 'J (z — fj,)(fi — X) J (z-X)(z-fi) 
Hence equation (I2.3P can be written in the form 



^Wfi( A ) dA = l [ p ™ ix) d \ i f &M^ (24) 



A f3n J (z — X) 2 n J (z — X)(z — fi) 

Let us introduce notations: 

k^'^\x, fj,)dXd/jL 



(z- X)(z- fi) 



(2.5) 



where 



fc&'W) = n(n - 1)P^{\H) - n 2 p^l(X) P ^%) + n8(X - /^(A). (2.6) 
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Moreover, we denote 

,fi,h 

A — z ' 1 ' z — X 

Then equation (|2.3j) takes the form 



1 /^(A), 1/2 \ /-/«(A) ;x 1 



n 



y s-a dA - » u " J y f^af " (2 - 8) 



Using that V(<2, C) is an analytic functions of £ in D, by the Cauchy theorem we obtain that 



if £ is encircling cr e , and z is outside of £. Thus (|2.8p takes the form 

i ^'(a)pK(a) 1/2 \ p£g(A) 1 



1/2 \ [ P\ fi ' h W i 

y ,-a dA - „ u " J y f^af ■ (2 - 9) 



On the other hand, it is easy to show that g„ satisfies the equation 

g 2 v (z) + V'(z)g v (z) + Q v (z) = 0, Q(z) = £ V(z, C)^(CK- (2-10) 

Hence 



= + ^ 2 U'(z) 2 -4Q„(z). 



2 v 7 2 

Using the inverse Stieltjes transform and comparing with (jl.23p . we get that 

2g v (z) + rjV'(z) = P v (z)X^ 2 (z). (2.11) 

where X r/ (z) is defined by (|1.24p for q=l. 

Write h = g v + n~ x u n ^. Then, subtracting (|2.1U|) from (|2.9|) and multiplying the 
result by n, we get 

(2g v (z) + r,V'{z))u n>r) (z) - ^- £ U(z, C)WCR = (2-12) 

where 

= y dA - (J - 1 j ^(z) + i^(z) j (2.13) 

Using (f2~Tll . we obtain from (l2~T2l) 



P^X^K,^) + Q n (z) = F(z), Q n (z) = ~2~jri V{z > O^nM^C- (2-14) 



Then, choosing C sufficiently close to [a v , &J to have all zeros of -P(C) and z outside of C, we 
get 

JL £ (W^ffl + QM) - F(0) n(c) ( C ,_ c) = o. (2.15) 

Since by definition (|2.14p Q n (C) is analytic in D, and z and all zeros of P are outside of C, 
the Cauchy theorem yields 

1 I QnitW _ Q _ 



2iri T c P v (()(z - 

Moreover, since 

u n , v (z) = ~(J Pi!&( X ) dX ~ J PvWd^j + nO(2 



z ) = nO(z ), z — > oo, 
we have 

X l r) / 2 {z)u n ^{z) = nO{z- 1 ). (2.16) 

Then the Cauchy theorem yields 

,1/2 



i r x v (CK„(C)rfC _ Y i /2( x ^ 



Finally, we obtain from (|2.15p 



= Wt / P fnd dC rv (2 - 17) 

2™AV (z) Jc p vKQK z -Q 



Set 



= dist{z,a e }, (2.18) 
Then for any z : d(z) = d equation (|2.17p implies 

u n,v(z) = 1/2 + 172—^ p 77^77 s , 2.19 

Xy 2 (z)P„(z) 2vr^ 1/2 (z) /£' ^(0(C " ^) 

where the contour £' contains z but does not contain zeros of P(z). It will be convenient for 
us to take £ as far from a e , as it is possible. According to [2j, for any j3 and rj we have the 
bounds: 



\ x (n,v), , Cnlogn Cn^log 1 / 2 /} , Cn 1 / 2 log 1/2 n , 

l<W (*)l < -51^-. K,(*)l < — ^ — - K„(*)l < — ^ — > ( 

where C is an absolute constant. Set 

M n (d) = sup \u n>r) (z)\. 

z:d(z)>d 

By (|2.16p and the maximum principle, there exists a point z : d(z) = d such that 

M n (d) = \u n>v (z)\. 

Then, using (|2.19p . the definition of F (|2,13p . and (|2.20p . we obtain the inequality 

M„(d)<^M + ^, c 2 <c (i + w\u 



2.20) 
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with some Co, C\ depending only on P and C of (|2.20p . Solving the above quadratic inequal- 
ity, we get 

" M n {d) > \ (dnd+ yJC 2 n 2 d 2 -4CiC 2 nlogn/d 4 ); 



M n (d) < Udnd - ^C 2 n 2 d 2 - 4CiC 2 ralogra/d 4 

Since the first inequality contradicts to (|2.20p . we conclude that for d > re -1 / 6 log n the second 
inequality holds. Hence we get 



,„(*)! < C lognd- 5 (z)(l + ||/i , || oo ), d(z) >n- 1 / 6 log 



n. 



(2.21) 



Note that the bound gives us that if we consider 92(A) = 3? (A — z) 1 or 92(A) = Q(A — z) 1 
with d(z) > n -1 / 6 logn, then 



n 



< ^n|b (s) ||oo(l + Halloo), 



(2.22) 



where V (s) is the s-th derivative of 92 (now we have s = 4 in view of (|2.2ip ) and 

w n = sCq log n 

We are going to use the following lemma, which is an analog of Lemma 3.11 of [12] . 



Lemma 1 If \2. 22\) holds for any real h : ||/i'||oo < A (A > 1) and some 92 such that 
Hv^lloo > Hv'lloo, then there exists an absolute constant C* such that for any ||/i'||oo < A/2 



^ ) (A J /iMA)^)dAdM<C^(l + A) 2 ||^W||^ e 



(2.23) 



The lemma was proven in [12], but for reader's convenience we give its proof here. 

Proof of LemmaUl Without loss of generality assume that w n > 1. Using the method of 
1121. consider the function 



Z n (t) = E v ,j3,h S exp 



t 

2ti 



n . 
n i=l J 



W n = tU„(l + A)\\<p 



Ml 



where ^n,p,h {•} is defined in (jl.5p with V replaced by V V: h of (|2.ip . It is easy to see that 

f^^(A i )-E Afc+t¥ , /M?n {^(A i )}) > ) 1>0. (2.24) 



d 2 

^2 logZ n (t) = (2u' n )" 2 E /3 ^ + ^ /2?i;n 



vl=l 



Hence in view of (|2.22j) . 

log Z n {t) = log Z n (t) - log Z n (0) = J log Z n (r)dr < log Z n (t) 

= |t|(2S n )- 1 E 7? ^ ife+t¥ , /2l5n |^ ^(Ai) - y ^(A)p^(A)dA^) I 

(4t W 2 in W - P„(A)) dA < t € [-1, 1]. 

Z n (t) < el*l < 3, te [-1,1], 



\t\n 
2w n 



Thus 
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and for any t 6 C, \t\ < 1 

\Z n (t)\ < Z n (m) < 3. (2.25) 



Then, by the Cauchy theorem, we have 

1 r Z n (t')dt' 



Km 



2nf wl=1 {t'-tf 



, 1 

< 12, \t\ < -. 
- , 2 



and therefore for \t\ < 



\Z n (t)\ = \Z n (0)- f^z' n {t)dt\> 1 -. 
Hence log Z n (t) is analytic for \t\ < ^j, and using the above bounds we have 



log^(0) = — <b l ^^dt<C. 



dt 2 2m ./j t | = i/24 t 3 



Finally, in view of (|2.24p we get 

J k nAh (X, fx)(p(X)(p(fj,)dXdfi = E vAh i I J2(<p(\i) - E r)i/3ih {<^(Ai)}) J > < ±Cw\ 



si=l 



□ 

Let us come back to the proof of TheoremQ] Applying the lemma to ip^ (A) = $l(z — A) -1 
and ip^(X) = ^s(z — A) -1 with d(z) > n _1 / 6 logn and using (|2.22p for ||/i'||oo < 2, we obtain 
for such z (cf flZgUD ) 

\Sff(z)\ < C log 2 nd~ w (z) (2.26) 

Then, using this bound and fj2 . 2 1 f) in (12.17P and taking into account that n ^ 
for d(z) > n -1 / 6 log n, we obtain that for ||/i'||oo < 1 

(2//3-1) i <?;(ck , „,,_ 5 . 



Applying Lemma [T] once more, we get Sp^ (z) < Cd~ 10 (z). Using the bound in (|2.19j) we 
prove (OHJ) . 

To prove (jl.29p consider the Poisson kernel 

".<*> - ^pta*,- < 2 - 28 > 

It is easy to see that for any integrable ip 

it) \t\\ 1 c- /" Vii^dii 



vr J n - (X + iy)' 
Hence we can use (|2.27p to prove that for \y\ > n~ 1//6 logn 

11^*^111 <Cy- 10 , for u n ^(X) := n(p^\x) - p v (X)), (2.29) 
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where ||.||2 is the standard norm in L2OR). Then we use the formula (see [12]) valid for any 
u G L 2 (R) 

f'OO r 

/ e-yy 2s - 1 \\V y *u\\ 2 2 dy = T(2s) / (1 + 2|£|)- 2 *|u(£)| 2 ^. (2.30) 
Jo Jr 

The formula for s = 5, the Parseval equation for the Fourier integral, and the Schwarz 
inequality yield 

/ <p(\)u n)V (X)d\ = — I <p{Qun, v (g)d£ 



< g((iMi a+ n y wii a ) //-" 2 y/ 2 

- rV2(2s) Uo II^*^|| 2 ^ • 

To estimate the last integral here we split it into two parts |y| > n _1//6 logn and \y\ < 
n -1 ' 6 log re. For the first integral we use (|2.29j) and for the second - (I2.20|) , 

(ii) To prove (jl.30p we start from the relations 

a v - a = (77 - l)a* + 0((r? - l) 2 ), b r) -b = (r } - 1)6* + 0(fa - l) 2 ), 
P v (z) - P(z) = ( v - l)p(z) + 0(fa - l) 2 ), 

g v {z) - g{z) = (rj- l)m(z) + 0(fa - 1) 2 |^|" 3 ), (2.31) 

where a*, 6* are some constant, p(z) is some analytic in D function, and m{z) is analytic 
outside [a, b\. The first line here follows from the results of [10J. The second line follows from 
the first one, if we use the representation (jl.25p . The third line follows from the representation 

g r) {z)= l -(X]/\z)P v (z)-r ] V'{z)) 

and the first two lines. Then (11.30P follows from (ll,28p . combined with (I2.3ip . 

(iii) Consider the functions Vnj of the form 

Vr,, t (X) = tr,V(X) + (1 - t)V<[ \\), (2.32) 

where V^ (0) is defined in (fOTT) . Let Q^\it) := Q n ,p\V % t] be defined by (JT2J) with V replaced 

by V Vt t. Then, evidently, = Qn,pW^]i an d Qn,p(ty corresponds to (see (jl.3ip ). 

Hence 

^logQg(l) - ^logQS(O) = ±£ cft^log Q$(t) (2.33) 

where (A;i) is the first marginal density, corresponding to V^. Using (|1.8p . one can 
check that for the distribution (jl.ip with V replaced by Vt the equilibrium density pt has the 
form 

Pt (\)=t Pv (\) + (l-t)$\\), pW (A) = ^lW (2.34) 



n 
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with , d„ of (ll.3ip . Hence using ()1.28f) for the last integral in (12.330 , we get 
logQ n AvV} = logQ n ^[V^}-n 2 ^£[V^}+n^£[ V V] 

'"' ' 1 dtl ( V V(z)-Vi°\z))u n , v (z,t)dz, 



2 (27T*) J Jc 2d 
Changing the variables in the corresponding integrals, we have 

\ogQ n ,p[V^] = l og Q; /3 +(^ + n(l-/3/2))log^ ! 

2 71 8 2 6 2 ' 

Then (OTj) follows. □ 

3 Matrix kernels for orthogonal and symplectic ensembles 

The orthonormal system {^i }&Lo satisfies the recursion relations 

A^(A) = + 6i"Vf(A) + 4^&(A), (3-1) 

which define a semi- infinite Jacobi matrix j( n \ It is known (see, e.g. [16] ) that 

l4 n) |<C, |6^|<C, \n-k\<en. (3.2) 

By orthogonality and the spectral theorem, we see 

(D&\ k = nsign(j - k)V'(J^) jk => (3.3) 
(^k* = 0, |j - fc| > 2m, |(r>W) i|fc | < nC, \j - n\ < nc. 

We are going to use the formula for S n> p, obtained in [23] (see also [8]). In order to present 
this formula, introduce some more notations: 

^1 - — Krn-2m+l> t t J n -2m+2' • • • ' t Pn-l) ' 

and 

M rs := (£*(»), ($(")) T ), D rs := (($(">)', ($i n) ) T ), 1 < r, S < 2. 
Observe that = 1 together with (.D^ )jfc = f° r |j — k| > 2m implies 

T n = l-D 12 M 21 . 

Then have (see [8]) 

S nA (X,fi) = K„, 2 (A,/i) + $i(A) T L» 12 e$ 2 ( M ) - $i(A) T G%$a0u), 

G := £>i 2 M 22 (l - D 2 iM 12 )- 1 D 21 (3.4) 
S„/ 2 , 4 (A, /x) = K n , 2 (A, /i) + $ 2 (A) T J Di 2 e$i( Ai ) - $ 2 (A) T Ge$ 2 (^), 
G := -Z> 2 i(l - M 12 D 21 )- 1 M 11 D 12 , 

where K n ^ 2 is defined in (ll.lOp . 

Proof of Theorem^ Using (|3.4|) . it is straightforward to see that the first assertion of Theorem 
[2] follows from the following lemma. 
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Lemma 2 Given any smooth functions f, g and any fixed A > there exists a C > such 
that for all n > 2m and all j, k £ {n — A, . . . ,n + A} one has 



(f^ n ) )(A) 5 (A)^ n) )(A)dA < ^(ll/lloo + H/lloo) (iMloo + Ib'IU); (3.5) 



(it) 6(/^ (n) )(A)|< -^(11/1100 + II/'IU); (tit) logdet(T n )>C. 



Indeed, taking in (i) / = g = 1, we obtain that all entries of which are used in (|3.4p are 
bounded by Cra -1 , hence, having that | detT^ 1 < C we obtain that G and G have entries 
bounded by nC. This proves representation (|1.18|) . 

Proof of Lemma® Assertions (i) and (ii) of Lemma [2] will be derived from the asymptotics 
of the orthogonal polynomials in Section HI We now prove assertion (iii), using Theorem [U 

Note that without lost of generality we can assume that a C (— 1, 1) and v* = in (|1.8p . 

Similarly to the proof of Theorem [T] we choose e in (I1.26D small enough to have all zeros 
of P of (jl.25p outside of a £ and use the results of [T7] that if we replace in all definition (jl.ip 
- (jl.4p the integration in M. by integration with respect to o e , then the new \ [V] will differ 
from Q n ,/3[^] by the factor (1 + 0(e _nc )), where c > does not depend on n, but depends on 
e. Moreover, the new marginal densities will differ from (|1.4|) by an additive error 0(e~ nc ). 
Hence starting from now, we assume that this replacement is made and all integrals below 
are in a £ . 

Consider the "approximating" function H a (Hamiltonian) (cf (jl.ip ) 

Q 

F a (Ax...A n ) = -n^^CAiJ+EloglAi-Aj^l^fAiJl^.fAj)) -n 2 £*, 
yW(A) = ^W(A), (3.6) 

a=l 

where V^ a ^(A) is defined in (jl,35p . and S* is defined in (jl.36p . Then 

ff(Ai...A n ) = fl (Ai...A„) + AH(Ai...A n ), Ai, . . . , A n G cr e , (3.7) 

n 

Aif(Ai...A n ) = ^loglA.-A.I J] l <Ta , £ (A l )l ffQ , i£ (A,)-2n^y(A,) + n 2 S*, 
V(X) = J> ffa , e (A)/" log|A-/i|p(/i)^. 

a =l Jo-\cr a 



Qnl= / e^^- A »)dA!...dA„. (3.8) 



Set 

By the Jensen inequality, we have 

H (AH) B . )fj < \ogQ n ^[V] - logQ<J> < /? (AH) H> p , (3.9) 

where 



< • • • >^:= Q;i /(• • • )^ ( ^A, < • • • > HatP := (Q^)- 1 1 ( 
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Let us estimate the r.h.s. of (|3.9p for (3 = 2. 

q 



(AH) H)f) = n(n-l) J pg(A^) log |A - M | fl ^ a JX)l (7al ^)-2n 2 {V,p^)+n 2 L*. 

Here and below (., .) means the inner product in L2(cr e ). But using the definition of V and 
L* we can rewrite the r.h.s. above as 

(AH) H ^ = £ j l. a ,A^ al ^)^g\X- li\kf\\ii)dXd^ (3.10) 
+n 2 ^L[l CTQ , £ (pg-p),l^(^-,); 



where the kernel fcg is defined in (I2.6P with rj = 1 and /i = 0. The term with 5(X — fj) 
from (|2.6p gives zero contribution here, because in our integration domain |A — fi\ > 5 (see 
(jl.26p ). The last inequality in (|3.10p is obtained as follows. Since dist {cr a>£ , o" Q ', e } > 8, and 
a £ C [—1, 1], we can construct 6-periodic even function L QjCr /(A) with 12 derivatives such that 

L a>a '(X - n) = log|A - n\, A G ct Qj£ , fj, G a a > >£ , L a ,a'(X) = 0, |A| > 5/2. (3.11) 

Then, 

oo 

L a ,a'(X - fj) = c k e ik < X -M s , X G a a>e , Li G a a ,, E , \c k \ < Ck~ 12 . (3.12) 

k=~ oo 

Note that we need 12 derivatives to have for the Fourier coefficients of La,a' the above bound, 
which will be used later. Hence 

J J-O'a.e ( A ) 1 ^ a ',e M lo S I A ~ ( A , 

CO . 

= £ c fc yi^,(A)l^ i ,( M )e^V3 e ^/3 fc (n) (A ^ ) ^ dA4 = 0(1)i 

fc=— oo 

Here we used the well known bound for /3 = 2 (see e.g. |17| ) 

p(A)^Mfc< n) (A,/x)dAd/z = y |p(A) - ^(/i)| 2 ir2 2 (A, M )dAd/x < CH^A)^. (3.13) 
For each term of the second sum in (|1.28|) we have similarly 

L^Ml-P^'M^-P)) ( 3 - 14 ) 

oo 
fc=— oo 

This follows from the results of [5], according to which for any smooth <p 

<p(\)(p { $(\) - P (X))dX = Oin-'WW^. (3.15) 
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To estimate the l.h.s. of (|3.9p . we study first the structure of Q^\- Since H a does not 
contain terms log |Aj — Aj| with Aj G cr a , £ , Xj £ cr a i :£ with a / a' , it can be written as 



felH \-k q =n 

„ fci n 

= / IIWAi)- II l ffe>e (A i )^-( Al -" A »)dAi...dA n (3.16) 

J=l J=fcl + -+fcg-l + l 



a=l 



where the "effective potentials" {Va } q a= i are defined in (jl.35p . Take {k* a } q a=l of f j 1 . 33 j) and 

set 

k\\...k*\ 

K k := j^jQk-,p/Qk*;p- (3-17) 



It is evident that 



1 9 fciH hfc 9 =n 

Choose e sufficiently small to provide that \^* a n/k a — 1| < ejj"' for a = 1, . . . , q, if /c*| < en, 
where is chosen for <t Qj£ and (n/k a )Va a \ according to Theorem [TJ 

Let us prove that there exist n, ^-independent C*, c* > and a£l' such that 

ae- c »" 2 ' |jfe-jfe*|>eh! (3 ' 18) 

The proof of the first bound is based on Theorem [TJ We write 

log (Wfc^/WQk.,/?) = ^(Q^/Q^^-log^^/Q^^) (3.19) 

+ log(<!Q^/WQ^ 

where Q* ka . p is defined in (TQ2D . Then (TOTD yields 

log (Q^p/QUp) = ^f{ £ n/k a + k a I(k a ) + O(l), (3.20) 
== ^ j[ dt f{]T V -(0 ~ Vff ka (0)un/k a ((,t)d< + (l - 0/2) log(d n/fcQ /2), 

where £ n/fca = £[(n/£; a )vj a) ] is the equilibrium energy dL6]), u n/ka (.,t), vff ka , and d n/fea are 
defined in ([OTP . Using (TOO]) and (I23TD . we get for I(k a ) of (^201) 

feal(fe«) - k* a I(k* a ) = b^ a \k a - k* a ) + Oil) + 0((k a - k* a ) 2 /n), (3.21) 

where the concrete form of b^ a ' is not important for us. Moreover, formula (II. 32ft for n = 
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k a = k* a + I and n = k* a yields 

log (k* a \Ql a ^/k a \Ql^) - / log (VVr(/3/2)) 



= -{j{K + l? + ^(1-/3/2)) log 
+ (^) 2 + y(l -m) log^ + ^logr(/3(A ; ; +J )/2) 



= 1(1 - 0/2) ( log(/3^/2) + I) - l^fc* I - + 0(i 3 /A£). (3.22) 
Here we used the Stirling formula 

loglW* + j)/2) = (p(k% + j)/2 - 1/2) ]og(P(k* a +j)/2) - + j)/2 + 0{n~ x ) 
[P(k* + j)/2 - 1/2) log03A£/2) + + i)/2 - 1/2) A _ p(k* a + j)/ 2 + o(J/*«) 



and the representation 



log = bg + ± - + o(/ 3 /^: 



Moreover, 



2 fc* 2(fc*) 



3fe^/3 (3A:*)^/3 3plk a 3/3/ 



+ 
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The last relation and (^T9l) - (ET221) yield 



log [K\Q ka ,l3/k a \Qk*,l3 



2 °n/k 



(3.23) 



where 



0/2) logn + 0(l), 



,(«) = (1 - /3/2) ( log(/3A;* /2n) + I) + log (v^/TCS/2)) + &( a > 



with of (|3.2ip . To estimate the difference of the energies we introduce the densities 

ft/toM = ^U^Pn/kJX), P {a) W = P(A)1. Q (A). (3.24) 

Then 
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n 2 [L 



Pn/k a > Pn/k a 



n 2 [L 



~(a) 

Pn/k a 'Pn/k 



~(a) (a) 
Pn/k^P-P 



Hence taking 



n 2 £^ :=n 2 {L[p^\p^] - (y,p^) + 2L[p^\p - p<«>]) : 



we obtain 



k a ^n/k a n 



n 2 [L 



_ Ja) _ „(«) 



+ 2L 



+2L 



< n z L 



t} ka -P {a \p-P {a) 
^k a -P [a) ^k a -P (a) 



(3.25) 



n 2 L. 



n/k a ■ 
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Here we used that in view of (|1.8jl and our assumption v* = 0, we have 



v(X) = J log |A - p\p(p)dp - V(A) = 0, A G cr a , v(A) < 0, A cr 
^l(A)-P (Q) (A) = pi Q / ) fcQ (A)>0, 



which implies 



(^pi/L-P (a) ) = (-^nl-P (a) )<0- 



Note that for any function </? : suppi/? C [— 1, 1], / <p = <p* 



L 



ip,tp 



L 



|</| 2 log2< -|^| 2 log2, 



p — ip*ip*,cp — (p*ip* 

where ?/>*(A) = 7r _1 (l — A 2 ) -1 ' 2 ].^^!], and we used the well known properties of tp*: 

J log | A — p\ip*(n)dp = — log 2, A G [—1, 1], J ip*(fi)dfj, = 1. 
Since, by (l334lh 



(3.26) 



(3.27) 



(3.28) 



we get 

n 2 L n/ka < - log2\k a /n - p* a \ 2 n 2 . 

In view of f|3.23j) and (|3.25p . to obtain the estimate for the difference of energies, it suffices to 
obtain the bound for (ka) 2 £ n /k^ ~ n 

2 £( Q ) from below. But it follows from (EOT]) and (fL23l 

that 

fi% a W-P {a) W = 0(K\ 1 / 2 ), v(X) = 0(K\ 3 / 2 ), AgV Q , (3.29) 

where 

d*a ■= K/ n - Pa- 

Since \d* a \ < In' 1 , the l.h.s. of (^261) is 0(\d* a | 3 ) = 0(n- 3 ) for k a = k* a . Hence 

" = n2L n/k* + Oin- 1 ). 

In order to estimate L n / k * , observe that if we denote v n / k * (A) the l.h.s. of fjl .8|) for (n/A;*)V, 



then we have 



where 



V*5 (A) -v(X) = d* a V^ (A) + v* /k . , A G a Q n supp jo^, . 



u n/fc» = SU P V*s( A ) = v n/k * a (c a ) = 0(d* a ), 



and c a is a middle point of suppp^,, . Moreover, the above argument implies that the 
contribution of the integrals over the domain <t Qj£ \ (a a n supp/r^, ) is 0(|d* | 5//2 ). Thus 



n 



L n/K =0{l), and KS n/ka - (k* a ) 2 £ n/K < -log2(k a - k*J 2 + 0(1). (3.30) 



*n2 



Then we obtain the first bound of (|3.18p from ()3.23|) and the last inequality, if take the sum 
with respect to a = 1, . . . , q and take into account that ^2 a (k a — A;*) = 0. 
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To obtain the second bound of (13.181) . we use the inequality 



(3k 2 

logQfc Q ,/3 < —^£ n /k a + Ck a \ogk a 
with some universal C (see [2j). Using the inequality for k a and A;*, we get in view of (|3.30p 
log (Qkcp/Qk&f) < -^log2\k a -k* a \ 2 + 0(nlogn). 

Hence we obtain the second inequality of (|3.18p . Now we are ready to find a bound for the 
l.h.s. of (EH). We have 



(3.31) 



fciH yk„=n 



fclH hk q =n 



k a k a/ L[p[ k ;\p{y ) ]-2nY,k a (v(\),p { *f 



Kq, K a i 

Then, for \k — k*\ < en we write similarly to (|3.12p and (|3.14p 



L 



\PW k n 



(a) JV) _ JL n W 



P" ',Pi$ 



]=-oo 
oo 



< 



/ 2 



j=-oo 



+ \{Pn /ka -^ a \e^) 



+ 



< 0(n~ 2 ) + C(n/fc a - (/i*)" 1 ) 2 + (n/k a/ - (p^)- 1 )'. 

To get the bound 0(n~ 2 ) for the first two terms in the r.h.s. here, we used (jl.29p for 
<Pk = e tknX / 3 and the bound for Cj from (|3.12p . and for the last two terms we used (|3.12p . 
combined with the estimate 

(Pn/k a ~ ^«),e^ A /3) < C\n/k a - (^)- 1 |. 

The estimate follows from (jl.23p and ()2.3ip , since these relations mean that p v is differentiable 
with respect to n at n = 1, and {prj) 1 ^ £ ^i[o"a,e]- For \k — k*\ > en, in view of the second line 
of (|3.18p . it suffices to use that the r.h.s. of (|3.3ip is 0(n 2 ). Thus we get finally 

\(AH) Ha ^\<C. 

The bound, (OU]) . and {fity . yield 

Qn,lQn/2,4 X " 



{Pn/k a , 



k, 



A/3) 



-1\ 2 



det(T n 



> 



Q n>2 (n/2)!2«; " 1 = 1 \Q K , 2 (k*j2)\2K 



Qk^,lQk*/2,A 
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Hence it is enough to prove that each multiplier in the r.h.s. is bounded from below. Consider 
the functions Vt of (|2.32[) . Then, as it was mentioned in the proof of Theorem Q] (iii), the 
limiting equilibrium density pt has the form (12.341 ). which corresponds to Vt. Hence V n j^t 
satisfies conditions of Theorem Q] for any t £ [0,1]. Moreover, if we introduce the matrix 
T n (t) for the potential Vt by the same way, as above, then T n (0) corresponds to GOE or 
GSE. Consider the function 

L(t)=logdetr fc .(t). (3.32) 
To prove that \L(1)\ < C it is enough to prove that 

|L(0)| < C, \L'(t)\ <C, t G [0, 1]. (3.33) 

The first inequality here follows from the results of [22J . To prove the second inequality we 
use (jl.20j) for V replaced by Vt. Then we get 

L'(t) = (fc*) 2 ( J AV(\)p i *i / t 2 \\)d\ + J AV(X)p[ k ^ t (X)dX - 2 J AV(X)p^(X)dx) , 

AV(X) = (n/k* a )V^(X)-V^l a (X). 

Using (jl.28p . we obtain that the first and the second terms of (|1 . 28[) give zero contributions 
in L'(t), hence L'(t) = 0(1). Thus we have proved the second inequality in (13. 33ft and so 
assertion (iii) of Lemma [2l As it was mentioned above, assertions (i) and (iii) of Lemma [21 
combined with imply that F$ and F$ of (fTTBD are bounded uniformly in n. 

To prove the first line of (11.37P for j3 = 1, we use (]1.18p 
J ' pt}{X)f{X)dX = n- 1 J S n , 1 (X,X)f(X)dX 

2m— 1 

P ^(x)f(x)dx + F ji ] / A%W^ [ :l k (\)f(\)dx. 

j,fc=-(2m-l) 

Thus we obtain the first line of (|1.37j) from (|3.15p and (I3.5[) (i). For j3 = 4 the proof is the 
same. 

The second line of (jl.37p follows from the first one in view of Lemma [TJ 

To obtain (jl.38j) . observe that we proved already that the l.h.s. of (jl.38j) is more than the 
r.h.s. Hence we are left to prove the opposite inequality. To this aim we use the inequality 
(13. 9p and the relation (I3.10p . Then each term in the first sum of (13.10p can be estimated by 
the same way, as for f3 = 2 by using the second line of (jl.37p instead of (|3.13p . Each term in 
the second sum in (I3.10p can be estimated by the same way, as for f3 = 2, if we use the first 
line of (071) . □ 

Proof of Theorem Convergence of the diagonal entries of K n> i and K n j 2 ^ to the 
correspondent limiting expressions follows from the first assertion of Theorem [2j (ii) of ()3 . 5f) 
and convergence of K n ^ to K^. The same is valid for 12-entries of K n ^\ and K n M^- Moreover, 
since eS n ^(X,fi) = — eS n ^(/j,, A), one has 

rn en 
(eS n ^){X,p) = - I S n ,i(t,n)dt, {eS n/2A ){X, p) = - / S n/2A (t, p) dt, 
Jx Jx 

which implies convergence of 21-entries of K n) \ and if n /2 ) 4-D 
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4 Uniform bounds for integrals with e{fip^i ) 

Set 

S„ = n- 2 / 3+K , 0<k< 1/3, (4.1) 
q 

a ±S n = \J 0~a,±8„, Va,±8 n = [ E 2a-1 T S n , E 2a ± <5n]> (°"a,-<5„ C CT Q C 0"a!,+<5n)- 
a=l 

Then, according to [5], we have 

^{X) = i? (A)co S nvrF n (A) (l + O^- 1 )), Ae U , 
V&(A) = i?i(A)sinnvrF ri _ 1 (A) (1 + 0(0), A G (4.2) 

where Rq(X), R±(X) are some smooth functions, which may behave like \X-^ 4 (X)\ near each 
E a (see (fQ4"|) for the definition of X). 

F n (X) = F(X) + -m (A), F(A) = p(p)d/i, F n _i(A) = F(A) + -mi (A), (4.3) 
n Jx n 

with /) of ()1.23p and smooth mo, mi, such that their first derivatives are bounded by |Jf- 1 /2(A)|. 
It will be important for us that 

aW J Ro(A)i?i(A)cos(7r(m (A)-mi(A))) = 1, (4.4) 
where is defined in (|3.ip . The relation follows from the fact (see [5]) that 

vraWn- 1 ((^)(A))Vt ) 1 (A) - ^(A)(^i(A))') = p(X) + 0(n~ l ) 
For | A — E a | < 5 n we have 

4»>(A) = n 1 / 6 ^ ) ^(n 2 / 3 d> a ((-ir(A-F a )))(l + 0(|A- J E a |)) 

+n- 1 /6 jB g) A ,-' (n 2 / 3 $ a ((-l)«(A - (1 + 0(|A - £ a |)) + O^" 1 ), 

ri n A(A) = n 1 / f3 ^ ) ^(n 2 / 3 <I> a ((-ir(A- J E a )))(l + 0(|A- J E a |)) (4.5) 
+n- 1 /6 jB g)Ai' (n 2 / 3 d> Q ((-l) a (A - E a ))) (1 + 0(\X - E a \)) + O^" 1 ). 

Moreover, 

l^ n) (A)| + |ri n) 1 (A)| <e- ncdist3/2 ^>, Ag1\m, 

Functions <3? a in (|4.2p are analytic in some neighborhood of and such that $ Q (A) = a a x + 
0(x 2 ) with some positive a a . 

The proof of Lemma [2] is based on the proposition: 

Proposition 1 Under conditions of C1-C3 for any smooth function f we have uniformly 
in <j +Sn 

e(M n) )(A) = f(X)R (X) C -^^l a _ Sn+Xn (X) (4.6) 
+ y( f(E a )-^^SL^ ( n 2 / 3 <K Q ( (_1)«(A - E Q 

ti\ ; (-i) a ^(o) v n v 

+0(n- 5 / 6 )) l |A _ Ba |<,5 n + er n (A) + 0{n~ l ), 
:= / Ai(t)dt, 
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where |Xn(A)| < n l l 2 C is a piecewise constant function which is a constant in each o- a -8 n 
and each interval (E a — 5 n , E a + 5 n ), and the remainder r n (X) admits the bound 

\r n (\)\d\ < Cn- 1 / 2 ^ 4 . (4.7) 

a +s n 

Similar representation is valid for e(/'0^ 1 ) if we replace sin by cos, Ro by R\, F n by F n _i, 
andB^ byB^. 

Proof. Let A £ a a -s n - Then, integrating by parts in (|4.2h . we obtain 

f(^\^=f(m^ nl, ' (ll] 



- / smnF„ (// — d\i. (4.8) 

Moreover, by (|4.5p . we have for |A — E?2a| < <5« 

^M/CMm = »-^/(A) »(^ A -^»-/("^-<-*'» 

E 2a -5 n <5>(\- F a ) 

+0(n- 5 /6) = w -i/a B g" ) /( £ 2a) * (re2/3 ^; ( ( A } " ^ 2a)) + (n-V2 )const + (n- 5 /6). 

Similar relations are valid for integrals near E2 a -i- Taking into account that the integrals 
over M \ a + s n are of the order 0(e~ ncSn ) = 0(e~ n3K/2 ), we obtain (|4.6|) . writing e(fip^) as 
a sum of the above integrals and similar ones (with integration from A to E2 a )- Then, for 
A 6 o- a _s n 

1 fE2a-l+8 n 1 fE2q+5n 



2 7£i-<5 n 2 J_B 2Q -i+5 n 

and r n is the sum of the terms, which are under the integrals in the r.h.s. of (I4.8p . 
Hence we are left to prove the bound for r n . Using that 

<(A) = (2vr)- 1 P(A)X 1 /2 (A ) + n^m'(A), \R(X)\ < C|X-V4(A)|, 

we have 



< 



d f^)R(p) 



d/j, + 0(n 



-5/6> 



dfi nF'(fi) 

< C(l + H/'Hoo)^ 1 ^ 3 / 4 = C(l + ll/'IDn- 1 / 2 - 3 */ 4 . 



□ 

Proof of Lemma\M Using recursion relations (|3.ip . it is easy to get that for any \j\ < 2m 

^•(A) = /o,(A)^ n) (A) +/i,(A)vt ) 1 (A), 

where _/bj and /ij are polynomials of degree at most Note that since a^, and &[ are 
bounded uniformly in n for k — n = o(n), /oj and have coefficients, bounded uniformly 
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in n. Hence assertion (ii) follows from Proposition [TJ Moreover, it follows from the above 
argument that to prove assertion (i) it suffices to estimate 

h := ( 5 V&,e(M n) )), h := G^l.^Kl)). J 3 := (g^\e(f^)). (4.9) 
with differentiable f,g. It follows from Proposition [T] that 

2q 



0=1 



where 



t -1; /" f/n a\ d mD M sSinnF n (A)sinnF n _i(A) 
A,o = n k n f(X)g(X)R (X)R 1 (X) — — dX 

h,2a = n- 1 / 3 B[l a) Bil a) f(E 2a )g(E 2a ) [ 2a+n ^(n 2 / 3 ^ 2a (X-E : 

JEo.„-8 n V 



2a, 

Ai (n 2 / 3 <5> 2a (X-E 2a )) 

and Ii,2«-i is the integral similar to I\ t2a for the region |A — £^a-i| < It is easy to see 
that 

T p (q) p (a) f{E a )g{E a ) 

h ' a ~ B ^ B ^ 2n(^(0))2 (1 + 0(1))l 

Moreover, (|4.7p and (ii) of Lemma [2] yield 

rfi-« , ») = -(rfi).''»)< / |e(5ri_ ) 1 )l|rn|dA = 0(n- 1 - 3K / 4 ). 
Hence we are left to find the bound for I\ q, 

Ji,o = (2n) / f(X)g{X)Ro{X)R 1 (X) — — dA 

<ro \-i / #m /up f\\r> ,^ cosn(F»(A) +F n _i(A)) , „ 

+(2n) / f{X)g(X)Ro{X)R 1 {X) — — dX = I 10 + J 10 . 



By dHU) and we obtain 



_ l + o(l) / 2 /(A)ff(A) 

2n y_ 2 P(A)XV2 (A ) dA 

Integrating by parts, one can obtain easily that = 0(n~ 2 <5 n 3 ^ 2 ) = 0(n~ 1 ~ 3K / 2 ). 
The other two integrals from (|4.9p can be estimated similarly. 



Acknowledgements. The author thanks Prof. T. Kriecherbauer and Prof. B. Eynard 
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